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Abstract 

By identifying a family of corner cutting schemes as a dimension elevation pro- 
cess of Gelfond-Bezier curves, we give a Miintz type condition for the conver- 
gence of the generated control polygons to the underlying curve. The surprising 
emergence of the Miintz condition in the problem raises the question of a pos- 
sible connection between the density questions of nested Chebyshev spaces and 
the convergence of the corresponding dimension elevation algorithms. 
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1. Introduction 

Let n be a fixed positive integer and let < ri < r2 < ...r„ < r„+i... < 
r™ < ... be an infinite strictly increasing sequence of positive real numbers. 
Given a polygon (Pqj -Pi: ^ri) in > 1, we apply the following corner 

cutting scheme : For i = 0,1,..., n, we set P° = Pi and for j = 1,2,..., we 
construct iteratively new polyg ons (^Pq , , — T Pn+j) using the inductive rule 



and for z = l,...,n + j — 1 



l-^Pl:^ + [i-^)pr (2) 

Figure [T] shows the first iteration of the corner cutting scheme on a planar 
polygon {Pq, Pi, P2, P3) for n = 3 and positive real numbers < ri < r2 < < 
r4. In the case the real numbers are given by = i for every index i, then 
we recognize the degree elevation algorithm of Bezier curves, and in which it 
is well known that the control polygons of the elevated degree converges to the 
underlying Bezier curve. Consider, now, the case in which ri = i for i — 1, ...,n 
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Figure 1: The first iteration of the corner cutting scheme JTJ and ^ for the parameters 
n = 3, ri = 2, r2 = 5, rs = 7 and r4 = 14. 




Figure 2: The sequence of polygons generated by the corner cutting scheme Q and (|2]l and 
parameters n = 3, ri = 1, r2 = 2, ra = 3 and rj = 2j for j > 4. (left, four iterations of the 
scheme; right, 100 iterations of the scheme). The red curve is the Bezier curve associated with 
the control polygon {Po, Pi, P2, P3)- 



and ri — 2i for i > n. Figure [5] (left) shows the generated polygons from 
the scheme ([1]) and ([2]) from four iterations, while Figure [5] (right) shows the 
generated polygons from 100 iterations. The figure suggests the convergence of 
the generated polygons to the Bezier curve with control points {Po,Pi, ...Pn)- 
Consider, now, the case in which Vi ~ i for i = while — for 

i > n. Figure [3] (left) shows the generated polygons from four iterations, while 
Figure [3] (right) shows the obtained polygons after 100 iterations. It is clear 
from the figure that the limiting polygon does not converge to the Bezier curve 
with control points {Pq, Pi, Pn)- As we will exhibit in this work, the main 
difference between the example of Figure[2]and the one of Figure[3]is the fact that 
in the former we have X^i^i V^i — 00, while in the latter we have X^i^i < 
00. Wc will show, as a particular case of our main result, that if = i for 
i ~ 1, ...,n, and limg^oo = 00, then the limiting polygon generated from the 
corner cutting scheme (H]) and ^ converges to the Bezier curve with control 
points (Pqj P\, Pn) if and only if the real numbers satisfy X^i^i l/'^* = 
The emergence of the limiting polygon as a Bezier curve in the case ri = i 
for i — l,...,ri can be hinted to as follows : the linear space formed by the 
monomials with exponents the numbers ri = i for i — 1,...,7t, and extended 
by a constant is given by = span{l,t,t^ , ...,t"); which is the linear space of 
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Figure 3: The sequence of polygons generated by the corner cutting scheme 1^ and ([2jl and 
parameters n = 3, ri = 1, r2 = 2, ra = 3 and rj = for j > 4. (left, four iterations of the 
scheme; right, 100 iterations of the scheme). The red curve is the Bezier curve associated with 
the control polygon (Pq, Pi, P2, Pa)- 

polynomial of degree n. The space E has a special basis (the Bernstein basis) 
in which the notion of Bezier curve can be defined. For general real numbers 
ri,i — and imitating the previous construction, we obtain the Miintz 

space F = span{l,t'^^ ,t^^ , ...,t^"). The linear space F also possess a special 
basis (the Gelfond-Bernstein basis) first defined by Hirschman and Widder Q 
and extended by Gelfond which is in a certain sense a generalization of 
the Bernstein basis to the Miintz space F (in the case = i,i = 1, the 
Gelfond-Bernstein basis coincide with the Bernstein basis). Using the Gelfond- 
Bernstein basis, we can canonically define the notion of Gelfond-Bezier curve 
with control points (Pq, Pi, Pn)- Now, consider, for example, the limiting 
polygon of the corner cutting scheme ((T]) and ^ for the case n — 3 and in 
which ri = 2, r2 = 4, ra = 5 and — 2i for j > 3. Figure S] shows the generated 
polygons from 100 iterations and also shows the Gelfond-Bezier curve associated 
with the Miintz space F span{l,r\e^,r^) = span{l,t'^,t'^,t^) and control 
polygon {Pq, Pi, P2, Pa)- The figure suggests that the limiting polygon converges 
to the Gelfond-Bezier curve. Therefore, The main objective of this paper is to, 
effectively, prove the following 

Theorem 1. Let n be a fixed positive integer and let < ri < r2 < ---Tn < 
r„+i < ... < r„j < ... be an infinite strictly increasing sequence of positive real 
numbers such that lim<,_).oo ''s = 00. Then the limiting polygon generated from 
a polygon {Po,Pi, ...,P„) in K*,s > 1 using the corner cutting scheme (QP and 
(0) converges (pointwise and uniformly) to the Gelfond-Bezier curve associated 
with the Miintz space {l,t^^ ,t^^ , .-.,t^'^) and control polygon (Pq, Pi, Pn) if 
and only if the real numbers ri satisfy the condition 



Let us contrast our theorem with the celebrated original Miintz theorem on 
the density of Miintz spaces [lo| 

Theorem 2. (Miintz Theorem) Let < ri < r2 < < r„+i < ... < 

< ... be an infinite strictly increasing sequence of positive real numbers such 
that lims_j.oo ?'s — oo- The Miintz space span(\,V^ , ...,t^"^, ...) is a dense subset 
of C([0, 1]) (the linear space of continuous functions on [0, 1] endowed with the 




(3) 
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Figure 4: The sequence of polygons generated from 100 iterations of the corner cutting scheme 
ni and ((2]| and parameters n = 3, ri = 2, r2 = 4,r3 = 14 and rj = 2j + 10 for j > 4. The 
red curve is the Gelfond-Bezier curve associated with the Miintz space span{l, t^, t^, 4^4) and 
control polygon {Po , Pi , P2 , Ps) 



uniform norm ) if and only if 



— = 00 
^ n 



The emergence of the Miintz condition ([3]) in both of Theorem [T] and Theo- 
rem [2] is rather surprising and may suggest a deep relation between the problem 
of density in Miintz spaces and the convergence of corner cutting schemes. In 
fact, as we will show in section 2, the corner cutting scheme ([1]) and (0) can be 
interpreted as a dimension elevation algorithm of Gelfond-Bezier curves. There- 
fore, Theorem [T] can be restated as claiming that under the condition that the 
sequence < ri < r2... < r„ < ... satisfies lims_i.oo = 00, the density of the 
Miintz space span(l,t^^ , ...,t^"^ , ...) is equivalent to the convergence of the di- 
mension elevation algorithm of Gelfond-Bezier curves to the underlying curve. 
We can push this analogy even further as follows : It has been shown in [l| 
that the Gelfond-Bernstein bases are limit of the Chebyshev-Bernstein bases 
of Miintz spaces over an interval [a, 1] as a goes to zero. From this property, 
it is not hard to show that the conditions of Theorem [T] are sufficient for the 
convergence of the dimension elevation algorithm of a Chebyshev-Bezier curve 
in Miintz spaces to the underlying curve. As the Chebyshev-Bernstein bases 
over an interval [a, b] can be defined for any linear space E = span{l, ui, Um) 
such that the space DE = span(u\ , ...,u',„) is an extended Chebyshev space of 
order m over the interval [a, b] ll|, l8[, we can ask for the following more general 
question : Let n be a fixed positive integer and let iti, W2, Wn, Um, ■■■ be an 
infinite sequence of C°° functions over an interval [a, b] such that for every fc > 1, 
the space Ek = span{l,ui,U2, ■■■,Uk) is such that DEk = span{u[,u'2.-.,u'i^) is 
an extended Chebyshev space of order k over the interval [a,b]. For any func- 
tion F 6 span(l, ui, U2, it„) with control polygon (Pq, P„) over the interval 
[a, b] , we can define the control polygons of the dimension elevation algorithm 
[9| with respect to the nested spaces En C C .... C E^ C the question 
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Figure 5: The sequence of polygons generated from 100 iterations of the corner cutting scheme 
(T) and lO and parameters n = 3,ri = l,r2 = 2,rs = 3 and for the left figure we have 
= 5 — i for j > i and for the right figure, rj = 50 — i for j > 4. The blue curve is the 
Bezier curve associated with the control polygon (Pq, Pi, P2, Ps)- 

is then 

Is there a connection between the density of the space Eoo = span{l, ui, 

U2, u„, Mm, ...)as a subset of C([a,b]) and the convergence of the (Q) 

associated dimension elevation algorithm to the underlying curve ? 

A hypothesis of equivalence is ruled out by the following fact: the condition 
lims_>oo = 00 can be dropped in Miintz theorem [5J however, such condition 
is necessary in Theorem [TJ For example, Figure [5] (left) shows the limiting 
polygon for the case n = 3, ri = 1, r2 = 2, = 3 and = 5 — i for j > 3. The 
limiting polygon does not converge to the Bezier curve with control polygon 
{Pq, Pi, P2, P3). As it will be clear within this work, the main reason for the 
non-convergence of the dimension elevation algorithm to the underlying curve in 
this case is the fact that the set of control points (?/™)o<i<m of the function f^^ 
with respect to the Miintz space span{l,t^^ ,t^^ , does not form a dense 

subset of the interval [0, 1] as m goes to infinity. It is interesting to note that for 
example when n — 3 and the real number are given by ri = 1, r2 = 2, = 3 
and ri = 50 — i , the limiting polygon is very close to the underlying curve and 
yet does not converge to the curve, as shown in Figure [5] (right). 

Regarding question (Q), we conjecture the following scenario : If for any 
fixed positive integer n, the dimension elevation algorithm with respect to the 
nested spaces En C i?n+i C .... C Em C ... over an interval [a, 6] converges to 
the underlying Chebyshev-Bezier curve then the space E^o is dense as a subset 
ofCi[a,b]). 

The proof of Theorem [T] consists in first showing, in section 2, that the 
corner cutting scheme ([1]) and Q can be interpreted as a dimension elevation 
algorithm of the Gelfond-Bezier curves. This allows us, in section 3, through 
a refinement of the elegant method of Prautzsch and Kobbelt [l^l to prove the 
theorem by induction on the fixed integer n. 

2. Gelfond-Bezier curves 

Let / be a smooth real function defined on an interval /. For any real 
numbers xq < xi < ... < Xn in the interval /, the divided difference [a;o, a^n]/ 
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of the function / supported at the point Xi,i — 0, n is recursively defined by 
N]/ ^ f{xo) and 



[xo,Xi,...,Xn]f 



[xi,...,Xn]f - [a:o,a:i,...,a;„_i]/ 



if n > 0. 



(4) 



If some of the Xi coincide, then the divided difference [xq, ...,Xn]f is defined as 
the hmit of (|4]) when the distance of the Xi become arbitrary small. A simple 
inductive argument show that when the Xi are pairwise distinct then we have 



i=0 ^j=0,j^i 
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f{xo) 
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Xl 
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fix^) 






V{xo, Xl 


Xn) 





(5) 



where V{xo, ...,Xn) is the Vandermonde determinant. Note that by the 
divided difference [ccq, xi, is symmetric in the arguments xq, xi, 
Consider, now, the function ft{x) — , where t is viewed as a parameter. For 
a sequence A = (0 — rg, ri, r„) of strictly increasing real numbers, we denote 
by i?A the Miintz space E\ = span{t^° ,t^^ , ...,t^"). 

Definition 1. For a sequence A — (0 = rQ,ri, ...,r„) of strictly increasing real 
numbers, the G elf ond- Bernstein basis of the Miintz space E\ with respect to the 
interval [0, 1] is defined by 



and 



H'^Ait) = i-^y rk+i...r4rk....,r„]ft for fc = 0,...,n-l 



The determinant representation of the divided differences ([5]) shows that for 
k = 0, ...,n — 1, the Gelfond-Bernstein basis can be expressed as 



rk+irk+2---rn 
y(rfc,rfc+i,...,r„) 





1 






1 


Tk+l 




1 





n-fc-1 
' k 

n-fe-1 
'^k+1 

„n—k— 1 



(6) 



Formula ([6]) reiterate the fact that every function H"^ k = 0,...,n is an ele- 
ment of the space E\. The Gelfond-Bernstein basis possesses several properties 
that are similar to the classical Bernstein basis over the interval [0, 1]. For the 
sequence A — (0, 1, 2, n), the Gelfond-Bernstein basis coincide with the clas- 
sical Bernstein basis. Moreover, for any sequence A = (0 = ro,ri^..,rn) of 
strictly increasing real numbers, and for any k = 0, n, we have |7| 



< Hl\{t) < 1 for any t G [0, 1] 
and for any t G [0, 1], we have 



(7) 



(8) 



k=0 



6 



Pi P2 




Po P3 



Figure 6: Gelfond-Bc5zier curves associated with the control polygon {Po, Pi, P2! P'i) and 
Miintz spaces : blue curve span(l,t,t'^ ,t^), red curve span(l, i, , i^"), green curve 

Moreover, the Gelfond-Bernstein basis is totally positive on [0,1], i.e. for any 
sequence < to < ti < ... < tn < 1, the matrix {H^ j^[tj))o<k,j<n is totally 
positive (i.e. all its minors are nonnegative). This property gives rise to the 
so-called variation diminishing property of Gelfond-Bezier curve, i.e. given a 
Gelfond-Bezier curve F with parametrization 

n 

P{t) = ^ HJ^^^{t)Pi with P, e M^ s > 1, (9) 

fe=0 

the number of intersections of any hyperplane in W with T does not exceed the 
number of intersection of the hyperplane with the control polygon [Pq, Pi, Pn). 
Note also that for the Gelfond-Bezier curve in we have P(0) = Pq and 
P(l) = Pn. Figure |6] shows examples of Gelfond-Bezier curves associated with 
a single control polygon {Pq , Pi , P2 , P3) and different sequences A. For a more 
thorough study of Gelfond-Bezier curves, we refer to [H. 

Lemma 1. Let Ai = (0 = tq, ri, r„) and A2 ^ (0 = tq, n, r„, r„+i) be 
two sequences of strictly increasing real numbers. Then, for fc = 0, ...,n 

HkM (t) = ""f "' gfe,!^ it) + 'i^Hltl.K. (*)• (10) 

Proof. From the definition of the Gelfond-Bernstein basis, the right hand side 
of equation (ITUl) is given by ( for fc < n — 1) 

(-l)""*-rfcrfc+i...r„ {[ru+i, ...,r„+i]/t - (r„+i - rk)[rk, ...,r„+i]/t) . (11) 

From the definition of the divided difference, we have 

[rfc+i, ...,r„+i]/t - [rfe, ...,r„]/t = (r„+i - rk)[rk, ■■■,rn+i]ft- 

Inserting the last equation into (fTT|) conclude the proof of the lemma for fc < 
n — 1. For k = n, the left hand side of (fTO)) is equal to 

- {rn+i - r„)[r„,r„+i]/t = t^- = Hlj,^{t). 

□ 
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Let Ai and A2 be the two sequences given in Lemma [TJ and let P be an 
element of the Miintz space i?Aj . As i?Aj C Ea^ , the function P can be expressed 
in both of the Gelfond-Bernstein bases associated with the two spaces as 

n n+1 

p{t) = ^iJ?,A, ^Y.^kXK^^P'^- (12) 

fc=0 fe=0 

Using Lemma [1] to detect the coefficients of H^'^{t) in the expansion (fT2|) . we 
readily find 

Corollary 1. The G elf ond-Bezier points Pk in are related to the Gelfond- 
Bezier points Pk by the relations 

Pq — -Poi Pn+l — Pn, (13) 

and for k — 1,2, n 

Pk = Pk-i +(l- Pk- (14) 

Note that equations and ([T^ describe the first iteration of the corner 
cutting scheme ([TJ and ©. Therefore, the corner cutting scheme can be can 
interpreted as an iterative dimension elevation of the Gelfond-Bezier curve P 
with respect to the nested Miintz spaces Ea^ C Ea^^i C •■•-Ea„ C where A^ 
refers to the sequence A^ = (0 = tq, ri, r^). 

For later use, we will need the following two lemmas, in which we omit the 
proofs as they can be readily obtained from the determinant representation ^ 
of the Gelfond-Bernstein bases. 

Lemma 2. Let H^^ k = 0, n be the Gelfond-Bernstein basis associated with 
the sequence A (0 = tq, ri, r2, r„). Then, for /c = 0, n we have 



rr"_ r 

j-rrn /.\ _ ilj—k+l 3 ijn+l f,\ 
ti^.,AW-|^^^^^(;-^i^.+l,A.W, 



where Ai = (0 = j'o, 1, ri + 1, r2 + 1, r„ + 1) 

Lemma 3. Let a be a positive real number, and let HJ} ,k = 0, (resp. 
HJ}f^^,k — 0, be the Gelfond-Bernstein basis associated with the sequence 

Ai = (0 = ro,ri,r2,...,r„) (resp. A2 = (0 = aro, ari, ar2, ar„). Then, for 
k = 0, n, we have 



3. The convergence of the dimension elevation algorithm 

The fundamental idea for the proof of Theorem [T] is essentially simple, and 
can be viewed as a refinement of the method of Prautzsch and Kobbelt [l^ . 
However, in practice, the simplicity of the idea is overshadowed by the com- 
plexity of the technical details. Therefore, to exhibit the fundamental idea of 
the proof, we will first use it for the proof the classical fact that the control 
polygons of the degree elevation of a Bezier curve converge to the underlying 
curve. 
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Let P be a polynomial of degree n represented in the Bernstein basis ,over 
the interval [0, 1], of degree m > n as 

m 

P{t) = Y.^TB^^\t). (15) 

By induction on n, we will prove that 

uv^\hT-P{-)\^0{-). (16) 
I m m 

For n < 1, we have 6™ = P{i/m), since the Bernstein-Bezier representation has 
linear precision. Now, let us assume that (|16p hold for polynomials of degree 
n—1. Given a polynomial P of degree n with the Bernstein representation ([T5]) . 
we consider the polynomial Q defined by 

Q{t) - P{t) - -P'(t). 

n 

The polynomial Q is of degree n — 1. Moreover, from ([T5|) . we have 

m , m — 1 

Q{t) ^Y^bTBrit) - - E cT-'Br-\t), (u) 

i=0 i=0 

where we have written the derivative P' as 

m— 1 



p'{t) = Y.cT-'Br-\t). 



i=0 

We can give explicit expressions for the coefficients c^~^ but, as we will see, 
such expressions will not be needed. Now, using the fact that 

we obtain from (fT7|) 

Qit) = b^B^it) + £ (C - i^CsA BTit). 
The induction hypothesis on Q shows that 

P{-)-bT)^ — (p'{-)-cTs^] -o(-). 

m J mn \ m J m 

The last equation leads to (using the fact that i < m) 

max\P{-)-bT\ = Oi-) + -m&x\P'{-)-ct~i^\. (18) 

Now, the polynomial P' is also of degree n—1 and therefore, we can apply the 
induction hypothesis on P', namely we have 

max|P'(^)-cr'|=0(- ^ ^ 



m — 1 m — 1 ' 
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We have 

Therefore, we have 

max |P'(-) - cTS^'l < I- - + 0{-) = 0(1), 

I 771 m m — 1 ni m 

where E = max[o^i] \P"{x)\. Insertmg the last equation into (fTS)) conclude the 
proof of (HH). 



Notations: In order to apply the previous idea to the case of arbitrary 
Miintz spaces, we will first set some notations. We define the following difference 
operator A acting on sequences as follows : If = (0 = ro, ri, r2, r^) is a 
sequence of strictly increasing real numbers, then 

AA„ = (0 = ro,r2 - ri,r3 - ri, ...,r„ - ri). 

For l<fc<m— l,we define the sequence A'^Am iteratively using the equation 

A'^'A™ = A(A''-iA„) with A"A™ = A„ 

Therefore, we have 

A'^A.n = (0 = ro,rk+i - rk,rk+2 - rk, - rk). 

Now, for i = 0, ...,m — k, we denote by 'r]^''\Am) the «th control point of the 
function ^'''=+1^'"'' with respect to the Miintz space E^k/^^, namely, we have 

r^^^\A^) = 0; ,;W,(A™) = 1 
and for i = l,...,m — fc — 1 

^f^(A™)- ft (1^) 

We will adopt the convention that if i < 0, then ?7,p''(Ao) = and also write 
r]f\Am) simply as ??i(A„). We have 

Theorem 3. Let A„ = (0 = rg, ri, r„) be a sequence of strictly increasing 
real numbers and let Am = (0 = tq, ri, r„, r„i) be a longer sequence of 
strictly increasing numbers. Let P be an element of the Miintz space E\^ written 
in the Gelfond- Bernstein bases of E\^^ and E\^ as 

n rn 

4=0 4 = 

Then, there exist (n — 1) constants Ck depending only on the function P and 
the finite parameters ri, ...,r„, such that 

n-2 

P{V^i^„^)'^'^')-bT\ <5]Cfc 
fc=0 

for all i = 0, ...,m. We adopt the convention that X)fc=o ~ ^- 



(20) 
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Proof. We will proceed by induction on n. For n — 1 and as rji{h.,ri) are the 
control points of the function f^^ with respect to the Miintz space i?A„ over the 
interval [0, 1], we have P(77i(Am)^'''"0 — 6™ = and the conclusion follows. Let 
us assume the claim of the theorem is true for any element of a Miintz space of 
order less or equal to n — 1. Let P be an element of the space _Ba„ written in 
the Gelfond-Bernstein bases of i?A„ and as 

n m 

Pit) = Y^p^H-^jt) ^Y^^THrAjt). (21) 

i=0 1=0 

Let us denote by A„ and A™ the sequences A„ = (0 = ro, 1, r2/ri, r„/ri) 
and Am = (0 = ro, 1, r2/ri, r„/ri, rm/ri). To the function P in (l2lT) . we 
associate the function P in the space i?A„ defined as 

n 

P{t)^Y.P^Hl^St). (22) 

i=0 

As the corner cutting scheme ([T]) and ^ associated with a sequence S = 
(ri, ...,r/) is invariant by a multiplication of every elements of S by the same 
scalar, we necessarily have 

m 

p(t) = 5]5ri/a„w- 

i=0 

Consider the following function Q defined by 

Q[t) = P{t) - ^tP'it). 

rn 

It can be readily checked that the function Q is an element of the Miintz space 
of order n — 1, Eyi„_-^ where n„_i = (0 = rg, 1, r2/r'i, r„_i/ri). Therefore, 
we can apply the induction hypothesis to the function Q. Before we apply such 
induction, let us first express the function Q in the Gelfond-Bernstein basis 
ff™. . We have 

m m — 1 

Q{t) = bTKA,„ it) - E ^r^'^ir.-. (^)' (23) 

2=0 " i=0 

where we have denoted P' as 

m— 1 



^'w = E^r-'^rr~-.W' (24) 

i=0 

where H™j7^ ^ is the Gelfond-Bernstein basis with respect to the sequence 
Tm-i = (0 = r-o,(r2/ri) - 1, (rs/ri) - 1, (r,„/ri) - 1). From Lemma [2j 
we have 

tiiixL.i^)= n i-^M^"i.Aji) = ^>+i(A™)i?r:i,Ajo- 
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Therefore, from 



we have 



i=l 



The last equation shows that we have apphed upon Q a dimension elevation 
from the Miintz space associated with the sequence n„_i to the Miintz space 
associated with the sequence Am- Moreover, it can easily be checked that 
Vi''\Am) = vif\h.m) for any i and A:. Therefore, the induction hypothesis and 
the expression of Q show that there exist (n — 2) constant Ct depending only 
on the polynomial Q and the parameters ri, r„_i such that for i — 0, m 
we have 



{P{V^{K^)) - hT) - ^miArn) (A,„) ) ) - C^') 



< 



n— 3 
fe=0 

Therefore, we have 

3 

|F(,7,(A„))-6r| <^Cfe 



^fl(A™)"(ry,|':|^|i)(A„0 



fc=0 

n 



'7a(A™)-(e|,^|i)(A„) 



(25) 



Now the function P' is an element of the Miintz space of order n — 1, £'r„_i, 
where r„_i = (0 = ro, (?'2/ri) - 1, {r^/ri) - 1, (r„/ri) - 1). The expression 
([M)) . shows that we have applied upon P' a dimension elevation from the Miintz 
space associated with the sequence r„_i to the Miintz space associated with 
the sequence Tm-i- Therefore, again by the induction hypothesis, there exist 
{n — 2) constants Ek depending only on P' and the parameters ri,...,r„ such 
that 



n — 3 



P'MVm-l)'^) 



k=0 



It can be easily shown that Ty^'^-' (r„i-i) = Vi'^^^^Am)- Therefore, we have 



71 — 3 



p'(^r^(A™)^)-crH<E^fc 



fc=0 



^fr^(A,n)-(ci.%i)(A™) 



(fe+2) 



Moreover, we have 

IP'MArn)) 



W-1 



< 



P'{V^{^m))-P'{{vt\i^r.))^) 

P'((ryli\(A„))^)-a 



W-1 



Therefore, 



\P'MA„,)) - < r7,(A„) - i4l\{A„.))^) 



C+ 



n-2 

E^- 

k=l 
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where C = max[Q i]\P" {x)\. Inserting the last inequahty into (1^ and using 
the obvious fact that |(?7i(A,„))| < 1, show that there exist (n — 2) constants Li 
depending only on the polynomial P and the real values ri , . . . , r„ such that 



n-2 



\P{miAm))-bT\<J2^k 



(26) 



fc=0 

In view of Lemma |3] and equation (|22p , we have 

Inserting the last equation into (pS)) conclude the proof of the theorem. □ 

The following lemma is implicit in '7], and even more explicit in [3], as our 
hypothesis are different from the ones taken in the latter and for the seek of 
completeness, we will include it proof. 

Lemma 4. Let ^ be a strictly positive and let aj and bj, j — 1, 2... be sequences 
of real numbers in ]0, 1[ such that 

li-?^-7. (27) 



j-i-oo In a 



■J 

us assume 



Define Ai{m) = Y['jLi+i % ^'^'^ Bi{m) — Hjli+i (* < '^^'^ 
that for any fixed i, we have 

lim Ai{m) = and lim B,(m) = 0. (28) 

Then 

lim {Ai{m)^ — Bi{m)) ~ uniformly in i. 

Proof. We should prove that for every ei > 0, there exist an mo such that for 
all m > mo we have 

\Ai{my — Bi{m)\ < ei for i = 1, 2, m. 

Let us fix an ei > and select an e g]0, 1[ such that 

1 — < ei, < ei; e < 7 and e^'"^ < ei. 

Condition (|27p shows that there exists a jo such that for any j > jo, we have 

7-e<^^^<7 + e. (29) 
log Oj 

Since logaj < 0, imply that for any j > jo, we have 

a] < bj < a] . 

The last equation shows that for any j > jo and for any m > jo we have 
Aj(m)T+'= < Bj{m) < Aj{m)'^-'. 
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We can rephrase the last assertion as follow : There exists an mo = jo such that 
for any m > mo, we have 

A, (mr+' <B, (m) < A,(mp-^ for j=joJa + l,...,m. (30) 

Let us fix m > mo = jo- If for a certain index j > jo, we have Aj{m)'^ > Bj{ni) 
and Aj{ni) < e, then, we have 

< Aj{my ~ Bj{m) < < ei. 

If for a certain index j > jo, we have Aj{m)'^ > Bj{m) and Ajlm) > e, then 
from and using the fact that Aj{m) < 1, we have 

< Aj{my-Bj{m) < Aj{my -Aj{m)''+'' = Aj{my {1-Aj{my) < < ex. 

If for a certain index j > jo, we have Aj{m)"' < Bjim) and Aj{m) < e, then, 
from we have 

< Bj{m) - Aj{my < Aj{my-' - Aj{my < e^"'^ < ei. 

Finally, if for a certain index j > jo, we have Aj{m)^ < Bj{m) and Aj{m) > e, 
then, from (PH)) . we have 

< Bj{m) - Aj{my < 1 - A] < 1 - < ei. 

As we have exhausted all the possible cases on the behavior of a pair of numbers 
Aj{m) and Bj{m) for a certain index j > jo, the conclusion of theses cases show 
that for any m > mo = jo , we have 

l^jirf^y - Bj{m)\ < ei for j ^ jo, jo + 1, ■■■,m. (31) 

Now condition (|28p. shows in particular that for any j < jo, there exists an 
Mo(j) such that for any m > Mo{j), we have 

\Aj{my - Bj{m)\ < ei. 

As we have a finite set of Mo{j),j = l,...,jo — 1, if we denote by Lo = 
maxj=i^... jQ„i(Mo(j)), then for any m > Lq, we have 

\A,(my ^ B,{m)\ < ex for j = 1, 2, jo - 1. (32) 

Therefore, by taking Mo — max{mo,Lo) and in view of ((3T|) and (|32|) . we have 
for any m > Mo 

\Aj{my — Bj{m)\ < ei for j = 1,2, ...,m. 

□ 

From the last lemma, we can prove the following 

Theorem 4. Let Aoc (0 — 7*0, T"!, r„, ■■■) be an infinite sequence of 
strictly increasing real numbers such that 

oo ^ 

lim Tg — oo and — = oo. (33) 

i— 1 ^ 
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For any integer m, we denote by A,„ the subsequence of A^o given by A„j = 
(0 — ro, ri, Tm)- Let P be an element of the Miintz space i?A„ written in the 
G elf ond- Bernstein bases of Ea^ and Ej^^ (m > n) as 

n m 

Pit) = ^p.i?;:A„(t) = Y.bTKAjt)- 

i=0 i=0 

Then 

lim P(?7,(A„0i/'^i) - 5™ = uniformly in i. (34) 

m— >oo 

Proof. In view of Theorem [3l we need only to show that under the conditions 
, the right hand side of ([20]) converges to zero as m goes to infinity, uniformly 
in i. As we have finite terms in the sum in (I20p . we will only need to show that 
for any fixed k such that < fc < n — 2, we have 

Jim^,7(!)jA,„)- (^(^+i|^j(A„0)'^'--+^"'^'=+^ -0 uniformly in (35) 



Let us first deal with the indices i such that z > fc + 1, in this case, if we denote 

by 

/, ?'/c+2 - ?'fc+i \ J , A rk+1 - rk\ . . „ 
a,- = 1 — and 6, = 1 ; ? > fc + 2, 

then, according to fjl9p . and imitating the notations of the Lemma [4l we have 

rn rn 

vt\i^o)= n b,=B,{m) and r;|^|,%(A„) = [] a, = A,{m). 

The fact that the sequence (rfc)o<fc<oo is strictly increasing, shows that aj and 
bj are elements of the interval ]0, 1[ (j > fc + 2). Moreover, as limj_).oo rj — oo 
and using the I'Hospital's rule, show that 

r Inbj rk+i - 

lim 1 = — -f > 0. 

j-+oolnaj rk+2 - rk+i 

To prove that for any fixed i > fc + 1, we have lim„i^oo ^i("i) — 0, we proceed 
as follows : Since aj s]0, 1[, we have 

l^^^rj^±l^JJ^ ^^dthcn -r^> ft (36) 
aj rj-rk+i Ai{m) V - rk+i J 

Moreover, 

= {rk+2 - Tk+i) V > (rfc+2 - rfe+i) V — . 

j=i+l 1 + i='i+l ^ ^ 

Therefore, the conditions (p3)) show that 

m m / \ 

lim _fe+2 k+i_ _ ^ thus lim TT ( 1 H — ) = cx). 
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which by conclude that linim^oo Ai{m) = 0. Similar treatments for bj show 
that for any fixed i > fc + 1, limm_^.oo Bi{m) = 0. Therefore, applying Lemma 
m (after an obvious shift of indices) shows the convergence of (IM)) uniformly in 
i > fc + 1. For i < fc + 1, then the term in is zero and if i = fc + 1 the 
term in (I35p is ri[''\Am), which converges to zero as m goes to infinity. Then, 
using the trick of finitude as at the end of the proof of Lemma El conclude the 
proof. □ 

In order to conclude the proof of the main Theorem [1] using Theorem 31 we 
need to show that the point set Dm — {7yi(A„j)^/''i , j = 0, ...,m} form a dense 
subset of the interval [0, 1] as m goes to infinity. For this aim, we need the 
following result proven by Hirschman and Widder Q and Gelfond 

Theorem 5. (T3irschman- Widder @], Gelfond 

Let Kao = (0 = ro, ri, r„, r;, ...) he an infinite sequence of strictly increas- 
ing real numbers such that 

lim r. = oo and > — — oo. 



To every continuous function f in the interval [0, 1] , we associate the operator 
B^°° (/) defined as 

m 

Bi-{f){x) =Y.f{ri,{Kmf''-')H:XJx), 

where Am = (0 = tq, ri, rm). Then the sequences Bm^{f) is uniformly 
convergent with limit f as m goes to infinity. 

Using the last theorem, we can now prove the following 

Proposition 1. Let A^o = (0 = rp, ri, r„, r;, ...) he an infinite sequence 
of strictly increasing real numbers such that 



oo ^ 

lim = oo and 7 — = 00 

i—l 

Denote by Dm the point set Dm = {"HiiAmY^^^ ,i — 0, ...,m}. Then, as m goes 
to infinity, the point set Dm form a dense subset of the interval [0, 1]. 

Proof. Let e be a strictly positive real number and let xq be a real number in 
the interval [0, 1]. Consider the continuous piecewise linear function / defined 
as 

-^ + 1 if 0<a;<2:o 

if xo < 2; < 1. 



fix) = 



\ \~xq \—xa 

As the function / is continuous, then by Theorem [SI there exist an mo such 
that for any m > mo, we have 

m 

I/(xo) - ^/(77.(A„Oi/^^)FX(^o)| < e. 
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Evaluating / in the last equation leads to 



2^0 



1 



(37) 



Now, if |a;o — ?7i(Am)^/''i | > e for i = l,2,...,m, then, using the fact that 
^,™g i?j"\^(a;o) — 1, shows that the left hand side expression of ([57| is also 
strictly greater than e, leading to a contradiction. Therefore, for any e > and 
any xo £ [0, 1], there exists an mo such that for any m > toq, there exists an 
i < m such that \xq — 7/i(A,„)-^/''i | < e. □ 

As this point, we are ready to prove Theorem [T] 

Proof of Theorem [1]: Let us prove the theorem when the sequence Aoo = (0 = 
To, ri, r„, r/, ...) of strictly increasing real numbers satisfy the conditions 
([32]). In this case, if we denote by P the Gelfond-Bezier curve with control points 
(Po, Pi, P„), we have to show that given a point x G [0,1] and a sequence 
of real numbers ?yi(a:) (Am)^^'^^ that converges to x as m goes to infinity (this 
is possible thanks to the density proposition [T]) , the points 6™^^ converges to 

P{x) as m goes to infinity. As the function P is continuous, P(77j(3.-)(Am)^/''i) 
converges to P{x). Therefore, for any e > 0, there exists an Mq, such that for 
any m > Mq, we have 

\P{x) ~ Pir^,^,^{A,^y/^^)\ < e. 

Moreover, from Theorem 01 there exists an Mi such that for any m > Mi 

|fe".)-^(^.(.)(A™)^/'^^)|<e. 

Therefore, \P{x) — ^ ^'^^ ^ — ^cix{Mq, Mi), thereby proving 

the pointwise convergence of the dimension elevated control polygons to the 
Gelfond-Bezier curve. The convergence is uniform as we have 

max |P(x)-6,^,) I < max |P(x)-P(r,,(,) (A„,)i/n ) |+ |p(^^^^^ (A™)^/'-^-^'!^.) I 

The function P is continuous in the compact interval [0, 1], thus 

max|P(a:;) — P(77j(j,)(Am)"'^^''^)| as m goes to infinity, 

and Theorem |4] shows that 

max |P(?yi(a;)(Am)^/''^) — b"L\ \ as m goes to infinity. 

This conclude the proof of the if part of the theorem. To prove the only if 
part of the theorem, we proceed by contradiction. Let us assume that the real 
number satisfy 

oo ^ 

— < oo. 

^ n 
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Without loss of generality, we can take the case in which the control points 
{Pq, Pi, Pn) are real numbers such that Pq = 0, Pi = 1 and Pi < P2 < •■• < 
P„. In this case, for any m, we have 6™ = and 6™ is a strictly decreasing 
function of m that converges to a strictly positive number < 5 < 1 



In this case we would have a gap between the point zero and 5, namely, if we 
take X = (5/2 than for any m > 0, we have 



and therefore, the limiting control polygon does not converge pointwise to the 
Gelfond-Bezier curve. 

4. Discussion 

In the following, we give a list of directions for future research as well as 
some open problems. 

1- The corner cutting scheme ((T]) and (HJ can be generalized so as to describe 
the dimension elevation algorithm of rational Gelfond-Bczier curves. In this 
case, weighted corner-cutting schemes are derived and the methods developed 
in this work can contribute to the study of the convergence of these new family 
of corner cutting schemes. 

2- We can study the limiting polygon of the corner cutting scheme in case we 
relax the hypothesis of strictly increasing sequences < ri < r2 < ... < r,„ < ... 
to only increasing sequences < ri < r2 <....< r„i < .... In this case, the 
Gelfond-Bezier curves involve logarithmic functions [5] , namely, if we rewrite the 
exponent (ri, r2, r„) = (fi, f2, r^) where the real number fi are distinct 
and if we denote by nij,j ~ l,...,m the number of indices i — 0, for 
which ri = ri, then the space spanll,^^ , ...,t^") = span{l, x"^^ (InxY; j — 
0,l,..m;i = 0,l,..., rrij — 1). The results of this work could be extended to this 
case by a limiting process. 

3- The Gelfond-Bezier curves are too "degenerate" at the origin to study 
the dimension elevation algorithm in case we impose no condition of mono- 
tonicity on the real numbers r^. For instance, if we consider the case n — 3, 
ri = 1, r2 = 2, ra =3 and rj — 1/ j, for j > 3 and we start with a control poly- 
gon (Pq, Pi, P2, P3) then the dimension elevated control polygon to the order m 
is not obtained by a corner cutting scheme similar to ([T|) and ^ but instead 
the algorithm collapses the first m — 3 control points to Pq while the remaining 
control points are given by (Pi,P2,P3) [ij. However, if we consider the dimen- 
sion elevation algorithm of Gelfond-Bezier curves far from the origin, i.e. over 
an interval [a, 1] with a > 0, then the Gelfond-Bernstein basis coincide with 
the Chebyshev-Bernstein basis the degeneracy at the origin disappear and 
the algorithm leads to a family of corner cutting schemes without imposing any 
condition of monotonicity on the real numbers r^. Unfortunately, such family 




\br-x\>- for 1 = 1,..., 



m 
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of corner cutting schemes involves rather comphcated coefficients expressed in 
term of Schur functions It will be interesting to find, for the far from the 
origin case, conditions on the real number for the convergence of the dimen- 
sion elevation algorithm to the underlying curve. In the theory of Miintz spaces 
over an interval [a, 1] with a > 0, and in which we impose no condition on the 
real numbers (beside that they are pairwise distinct), then the corresponding 
Miintz space is a dense subset of C([a, 1 
satisfy the so-called full Miintz condition 



if and only if the real numbers 



y- 



(38) 



The question is then does the surprising emergence of the Miintz condition in 
Theorem [T] for the real numbers with the condition of the theorem, repeat 
itself for the full Miintz condition (1551) for the far from the origin case. 

4- It is not difficult to show that with the conditions of Theorem [1] the con- 
ditions (j33|) are sufficient for the uniform convergence of the dimension elevation 
algorithm of Chebyshev-Bezier curve of the associated Miintz space to the un- 
derlying curve (the proof will be published elsewhere). However, the pointwise 
convergence is more involved and the question of whether the Miintz condition 
is necessary prove to be interesting. 

5- It is probably a difficult problem to study the rate of convergence of 
the corner cutting scheme ((T|) and ^ in case the real numbers satisfy the 
condition of Theorem [TJ Adapting the method of Prautzsch and Kobbelt 

to this problem shows for example that if the numbers are integers and that 
there exists a constant K such that Vj < K + j for all j > 1 then the rate of 
convergence of the corner cutting scheme is in 0{:^). 

6- It may happen that studying the limiting polygon of the corner cutting 
scheme (1) and (2) is richer under the condition 

CO ^ 

fe=i 

in analogy with the problem of studying the uniform closure of the Miintz space 
Eoo = span{l,t^^ , ...,t^" , ...) under the condition (p9|) . 
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